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I. INTRODUCTION 

Semiconductor solid state structures are believed to be 
the most promising prospective hardware for building a 
quantum computer Their decisive advantage, based 
on the existence of the huge semiconductor electronics 
technology, is their scalability*"^ once few-qubit opera- 
tions are successfully accomplished, the well-established 
modern semiconductor microelectronics industry greatly 
facilitates scaling up to many qubits. Although longco- 
herence times have been demonstrated in both GaAs^HIS] 
and Sl^tL^ semiconductor devices, the two architectures 
are at different stages of development en route to an 
eventually viable quantum computer. In the GaAs quan- 
tum dot system, one- and two-qub it m anipulations have 
been experimentally demonst rated EEHHli] However, in its 
counterpart silicon-based system, these manipulations 
are yet to be satisfactorily achieved, in spite of much 
efforts in several laboratories across the world. 

Several features distinguish the silicon-based devices 
from the GaAs systems. First, the silicon quantum 
dot systems generally have longer spin coherence time 
than the GaAs system, due to some novel properties 
such a s wea k spin-orbit coupling and low nuclear spin 
densityP^HMUl In this regard, silicon-based devices are 
often considered to be ideal for quantum computation, 
because quantum error corrections protocols are easier to 
implement for qubits with longer coherence times. How- 
ever, various traps and impurities introduced during the 
fabrication lead to unpredictable local fluctuations of the 
confinement potential, which has so far prevented a co- 
herent manipulation in these quantum dot systems In 
this regard, the Si quantum dot qubit is inferior to the 
GaAs quantum dot structure since the interface quality in 
the Si system, particularly for the Si/Si02 system (but 
also to a lesser degree for the Si/SiGe system), is sub- 
stantially worse than the corresponding epitaxially per- 
fect GaAs/AlGaAs system. Second, the valley degree of 
freedom resulting from the special band structure of sili- 



con, with bulk Si having six equivalent conduction band 
minima (i.e. a valley degeneracy of six), further compli- 
cates the problem. \y e also note that in silicon metal- 
oxide-semiconductor (Si-MOS) system as well as the SiGe 
system the Coulomb interactions are stronger than that 
of the GaAs devices, due to lower dielectric constant and 
larger effective mass, which consequently affects the ex- 
change interaction. 31 The larger electron effective mass 
of Si compared with GaAs makes the quantum tunnel- 
ing weaker and the effective Bohr radius smaller, leading 
to tighter lithographic constraints in Si over GaAs for 
achieving equivalent qubit properties, thus rendering the 
fabrication of Si qubits much more difficult experimen- 
tally. 

In spite of these experimental difficulties, considerable 
progress has recently been made experimentally in si li- 
con based systems, incl uding both Si-MOS devices^^HM] 
and Si/SiGe structures' 16 1 27 1 28 1 3 ^"^ Particularly impor- 
tant in the few-electron regime of a d ouble quantum dot 
system is the charge stability diagra m 39 1 40 1 (explained in 
more detail below), which is an intuitive and physical 
manifestation of the Coulomb blockade 41 physics as well 
as the controllability of the qubit via, for example, the 
tunnel coupling between the two dots. The variations 
of the tunnel coupling not only affect the exchange in- 
teraction in the actual qubit manipulation in a chosen 
subspace of the charge stability diagram, but also lead 
to a change in the geometry of the overall charge sta- 
bility diagram. Lai et a/P^ and Simmons et aZ™ have 
demonstrated highly tunable double quantum dot sys- 
tems on Si-MOS and Si/SiGe systems, respectively. The 
charge stability diagrams over a wide range of electron 
occupancy were obtained in these experiments, and most 
interestingly they have shown how the patterns of the 
stability diagrams change as the tunnel coupling between 
the two dots is varied. In the experiment on Si/SiGe^ 
the numerical value of the tunnel couplin g is extracted 
from experiments using a two-level model. 44 These ex- 
periments have paved the way for the realization of sin- 
glet/triplet qubits in silicon, and a quantitative under- 
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standing of quantum effects in the charge stability di- 
agrams in Si coupled quantum dot experiments is the 
central goal of the present paper. 

In a previous publication^ we have advocated a quan- 
tum theory using the generalized Hubbard model to de- 
scribe the double quantum dot system. The advantage of 
our theory over the traditionally used capacitance circuit 
modeP^ED ^ our moc [ e i naturally includes quantum 
fluctuations, which is particularly suitable for the discus- 
sion of the physics of the tunnel coupling in the context 
of charge stability diagrams. Although our approach is 
completely general since the only restrictions are sym- 
metry considerations, we focused on the GaAs system in 
our original work. In this paper we apply our method 
to the silicon quantum dot system. We start with a qual- 
itative survey of the experimental results, followed by a 
detailed explanation of the assumptions and simplifica- 
tions that we have made in applying our general theory 
to the Si double quantum dot system. Then we demon- 
strate that upon appropriately choosing the parameters, 
our theory is capable of giving quantitative description 
of the experimental charge stability diagram and tunnel 
coupling in the Si system. We also outline several future 
research directions implied by our work. 

The remainder of the paper is organized as follows. 
In Sec. [IT] we briefly discuss experimental results and 
assumptions that have been employed. In Sec. Ill we 
explain the generalized Hubbard model and the micro- 
scopic theory used in this work. Sec. [IV] presents cal- 
culated results on the charge stability diagram for both 
Si-MOS and Si/SiGe system, and the tunnel coupling in 
the Si/SiGe system. The results are compared to exper- 
iments side- by-side. Sec. |V| includes our discussion and 
conclusions. 



II. QUALITATIVE CONSIDERATIONS 

In the typical experimental setup of Ref. [42] the dou- 
ble quantum dots are confined by three barrier gates (left, 
central and right) and the depth of the potential wells in 
the dots is controlled by two plunger gates with voltages 
Vpi and Vp2 respectively, in analogy to the left (Vl) and 
right (Vr) gate voltages of GaAs devices. 18 To shrink 
the dot size, the plunger gates are eliminated in Ref. [43j 
and the gate voltages Vl and Vr are used instead. The 
Coulomb blockade^ means that the electron occupancy 
in each dot (denoted by rii for dot i, i = 1, 2) is typically 
changed by an integer as the gate voltages are swept, 
which can be directly detected by measuring the tun- 
neling current, using a quantum point contact.^ As a 
consequence, a plot of the differential conductance on 
the plane with Vp± and Vp2 (or Vl and Vr) as axes is 
a visualization of electron configuration in the quantum 
dot system and is therefore called the charge stability di- 
agram. Fig. [TJa)-(c) show the charge stability diagrams 
measured by Lai et at. on Si-MOS devices, reproduced 
from Fig. 2(a)-(c) of Ref. l42l The white regimes indicate 



the Coulomb blockade plateaus in which some particular 
set of electron occupancy (711,712) is maintained; the col- 
ored lines separating the white regimes represent the elec- 
tron tunneling changing the occupancy by integer num- 
bers. Fig. J2ja)-(c) show similar results in Si/SiGe struc- 
tures, reproduced from Fig. 2(c)- (a) of Ref. [43j [Note 
that the order in the two sets of figures is reversed here 
to comply with the convention that inter-dot coupling 
increases from panel (a) to panel (c).] 

The tunnel coupling between the two quantum dots 
can be adjusted by changing the electrostatic potential 
of the central barrier gate, which consequently affects the 
geometry of the charge stability diagram. The variation 
from panels (a) to (c) of Figs, [I] and [^illustrates t his pr o- 
cess, whose qualitative feature is well understood ! 46 * 47 * If 
the central barrier is high enough to well separate the two 
dots (so that the coupling between them is weak), the 
two dots behave basically independently and the charge 
stability diagram has checker board pattern [Fig. [TJa)]. 
As the central barrier is lowered the inter-dot coupling 
becomes intermediate, the charge stability diagram be- 
comes honeycomb-shaped [Fig. [TJb), Fig. [2|a) and (b)]. 
In Fig. [TJc) and Fig. [2jc) the inter-dot coupling is strong 
due to further lowering of the central barrier such that the 
two dots effectively merge into one, the stability diagram 
becomes like a series of diagonal parallel lines separating 
regimes with different total occupancy N = 711+712, while 
the boundaries between different (711,712) states within 
the same N block [e.g. (711 + 1,712) —> (711,712 + 1)] are 
hardly visible due to the hybridization from the quan- 
tum fluctuations. In this situation with inter-dot tun- 
neling induced strong quantum fluctuations, the classical 
assumption of integer electron filling per dot becomes 
invalid, thus making the classical circuit model inappli- 
cable. 

The prevailing quantitative approach of understand- 
ing the charge stability diagram is the capacitance cir- 
cuit model, which symbolizes the Coulomb interaction 
by capacitors. 46 In our previous worlP^ we have shown 
a mapping of the capacitance circuit model to the spe- 
cific limit of the generalized Hubbard model with quan- 
tum fluctuations completely suppressed. In Fig.[TJa), the 
inter-dot coupling including both Coulomb interaction 
and the quantum fluctuations is small, and the capaci- 
tance model is, to a good approximation, valid. However, 
in Fig. [ljb) the boundaries separating different states are 
rounded near the triple points due to quantum effects, 
while in Fig. [TJc) the quantum effects (rounding) become 
even more evident. The rounding feature of boundary 
lines cannot even be qualitatively explained in the ca- 
pacitance model in any way and the basic purpose of 
the present paper is to give a quantitative description 
of the experiments using our proposed theoretical ap- 
proach, combining the generalized Hubbard model and 
the microscopic theory. 

In principle, a detailed quantitative theoretical descrip- 
tion of a specific experiment needs to incorporate as many 
experimental parameters as possible extracted directly 
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from the data. However in our situation several com- 
plications arise. These Si quantum dot experiments are 
done in the multi-electron regime, and the exact electron 
occupancy (711,712) is unknown. While the capacitance 
model predicts a periodic form of the charge stability di- 
agram, the actual situation is far from clear (and does not 
seem periodic). A realistic multi-electron and multi-band 
calculation would be extremely time-consuming and com- 
pletely impractical. Therefore we shall focus on the 
"effective" two-electron regime, meaning that we select 
one diamond in each stability diagram [indicated as blue 
square frames in Fig. [TJa)-(c)], and argue that the elec- 
tron occupancies in the regime enclosed by the frame 
ranging from (711,712) to (ni + 2,ri2 + 2) are approx- 
imated by effective electron occupancies ranging from 
(0,0) to (2,2), although in reality the actual configura- 
tion is changing from (711,712) to (ni + 2, 712 + 2) with 
both ni and 712 being unknown (and possibly of order 
10-20). We shall further assume that the multi-electron 
multi-band effect can be safely regarded as a renormal- 
ization of the effective two-electron and one-band Hub- 
bard model parameters such as the Coulomb interaction 
U but otherwise the same physics remains. 49 Indeed a 
close examination of the experimental plots Fig. [TJa)-(c) 
shows that the stability diagrams are, to a good approx- 
imation, periodic over a wide range, with details varying 
slightly as the occupancies are changed. This indicates 
that our assumptions are reasonable. A complete clar- 
ification of this multi-electron issue requires further ex- 
perimental investigation on a much wider range of the 
stability diagram, new experimental techniques that en- 
able precise measurement of electron occupancies, and a 
substantially more complicated theoretical and computa- 
tional study on a multi-band system. This is far beyond 
the scope of the present paper and probably the current 
experiments. 

Another consequence of the aforementioned approxi- 
mation is that some of the gate voltages (for example 
the plunger gate voltages Vpi and Vp 2 ) are no longer ba- 
sic variables in our theory since they control all electrons 
while electrons in the effective two-electron regime sim- 
ply feel a different potential which depends on the overall 
electrostatic potential, the separation between the fabri- 
cated gate and the actual localized electron gas, and the 
underlying Fermi sea of the background two-dimensional 
electron gas in the semiconductor heterostructure. To fa- 
cilitate the discussion we adopt the language that is used 
in the study of GaAs devices, i.e. the left and right gate 
voltages Vl and Vr as the effective gate voltages for Si- 
MOS devices. We will not try to establish a precise con- 
nection between (Vp2, Vpi) and (Vr,Vl) partly because 
one has no reliable information of the electrons in the 
underlying Fermi sea and partly because the fabrication- 
induced disorder may come into play. We shall, however, 
assume that are related to (Vp 2: Vpi) linearly, 

i.e. the calculated stability diagram should fit the exper- 
imental plot after appropriate rescaling and shifting. In 
fact, we have slightly rescaled the ^-direction of the ex- 



perimental plot of Ref.l42lin Fig.JTJa)-(c) since according 
to our previous work 45 the boundaries separating differ- 
ent states within the same N block [for example the one 
separating (1, 1) and (0,2) regime] should be at 45° with 
respect to the x-axis (i.e. parallel to the line Vl = Vr) 
if the mapping to the capacitance model is valid. The 
experimental plots shown in Ref. [42] need be compressed 
along the ^-direction in order to exhibit this property. 
This rescaling is simply a matter of convenience and con- 
vention, and does not affect any of our results or conclu- 
sions. 

In addition to the assumptions above, we further as- 
sume that the Si conduction band valley splitting is suf- 
ficiently large at the interface so that the valley degree of 
freedom can be disregarded. The neglecting of the val- 
ley splitting because of interface effects is the standard 
approximation in essentially all of the Si spin qubit liter- 
ature with only a few exceptions ! 3Q | 5Q | We focus our study 
on the case without any external magnetic field, but the 
generalization to the finite magnetic field situation^ is 
straightforward. 

III. METHODS 

In this section we discuss the methods used to model 
and quantify the experiments. We first describe the gen- 
eralized Hubbard model, and then the microscopic theory 
where we use two different models of confinement poten- 
tials since the actual quantum dot confinement potential 
is unknown. We also briefly outline how we fix the pa- 
rameters to best describe the experiments. 

In our previous publication}- 4 ^ we have proposed a gen- 
eralized form of Hubbard model, retaining all terms al- 
lowed by symmetry (total particle number N and total 
spin S z conserved). The Hamiltonian can be expressed 
as 

H = H^ + H t + Hu + Hj. (1) 

The one-body part includes the chemical potential/level 
energy terms 

Hjjl = — HiTiifj, (2) 

and the tunnel coupling (hopping) terms 

H t = - J2 (^Lc2a + fc.c.) , (3) 

where rii a = c\ a Ci a , \±i is the chemical potential of elec- 
trons on site i (i = 1, 2), and t is the tunnel coupling or 
inter-site hopping. 

The two-body part consists of a direct Coulomb repul- 
sion term Hjj 

H v = Uin^n n + U 2 n 2 ^n 2i + ^12(^1^2; + n u n 2t ) 
+ (U 12 - J e )(ni t n 2t + n n n n ), 

(4) 



4 



and a term Hj including spin-exchange (J e ), pair- 
hopping (Jp), and occupation- modulated hopping terms 
(Jt) 

Hj = - j e c^4t C2 i Cl t _ J P c U c li c n c n 

ia 

The gate voltages Vl and Vr appear in the generalized 
Hubbard model i n term s of the chemical potentials /ii ? 2 
with the relatiorPEMl 



fJLi = \e\[a 1 V L + {l-a 1 )V R ] 
= \e\ [(1 - a 2 ) Vl + ot 2 V R ] 



-71; 

"72- 



(6) 

If the mapping to the capacitance modeP^ is valid, we 
have 



(U2 - U12) Ux 



TJ 2 



OL2 



(l/i - C/12) U 2 



UtU 2 



r/ 2 



(7) 



We note that when the quantum fluctuations are notice- 
able such that the capacitance model does not apply, one 
is not necessarily restricted to Eq. ^ to convert the 
result from (/i 2 ,/ii) to (Vr, V^) plane. For simplicity, 
in this work we shall only use the Coulomb parameters 
U\, U 2 , U12 an d the tunnel coupling t to fit the experi- 
ments, assuming other exchange parameters to be small, 
although they are invariably present in both microscopic 
calculations and actual experiments. Our motivation at 
this initial stage of the development of the subject is to 
keep as few free parameters as possible in the theory and 
still incorporate quantum fluctuation effects not accessi- 
ble in the widely used capacitance circuit models. 

As discussed in Ref. 45, the generalized Hu bbard model 
must be backed up by a microscopic theoryP^^Hl The 
many-electron Schrodinger equation is solved in the pres- 
ence of a double-well confinement potential. In each 
well the potential is approximated as a harmonic oscilla- 
tor whose low- lying eigenstates are Fock-Darwin states. 
According to the Hund-Mulliken approach, the two- 
electron wave functions are constructed by hybridizing 
and orthogonalizing the Fock-Darwin states. These two- 
electron wave functions span a Hilbert space from which 
all the parameters of the generalized Hubbard model are 
derived. The details are described in Ref. 45. Here, we 
use two different kinds of potentials: a biquadratic form 
[Eq. ^] and a Gaussian form [Eq. ([13])], and make a 
side-by-side comparis on of th e results. We start with the 
biquadratic potential ) 31 1 59 1 60 1 
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{TYIUJ 
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2 + y 2 ]-Hu 



-[{x-a) 2 + y 2 \ -M2,0}. 



(8) 



The biquadratic potential has a simple form which only 
involves three parameters: uo 2 and the inter-dot dis- 
tance 2a. We expect that the on-site Coulomb interac- 
tions U\ and U 2 are directly related to uj\ and uj 2i while 



the inter-dot Coulomb interaction should vary primarily 
with the inter-dot distance 2a with possible minor de- 
pendence on uji and uj 2 . This follows from the "local" 
assumption of the electron wave function that the elec- 
trons are well localized in their own respective potential 
wells. This assumption is built into the Hund-Mulliken 
approximation, where the many-body wave functions are 
constructed using local wave functions. (The constructed 
many-body wave functions are the true ones when the 
two potential wells are far from each other, or the central 
potential barrier is high, such that the two dots behave 
independently.) To understand this better we consider an 
isolated single quantum dot with harmonic confinement 
potential 



Vi (r) = irac^ 2 



(x - Xi) + (y - ViY 



(9) 



The subscript i = 1,2 labels the dots. If the lowest or- 
bital is occupied by two electrons with opposite spins, the 
Coulomb interaction between them (i.e. the two electrons 
in the same dot labeled by "i") is 

Ui = J dri J dr 2 \(fi (ri)| 2 |^ (r 2 )| 2 



■ r 2 \ 



= ke\r- 1 -= ke* J 7 ^ 
V 2 L V 2h 



(10) 



where k = 1/(4ty£q£) : the ^-orbital wave function is 



exp 



(x - xi) + (y - ViY 
21? 



(11) 



and the effective length Z, = y/h/ (mwj). Therefore, the 
on-site Coulomb interactions only depend on u>i in the 
isolated dot limit. If we ignore the overlap between Fock- 
Darwin states of different dots, the inter-site Coulomb 
interaction is 

U12 = J dri J ' dr 2 \<pi (ri)| 2 \p 2 (t"2)| 2 



r 2 \ 



= he' 



If + 11 



2a 2 



q + q 



2a 2 



q + ti 



(12) 



Here 2a is the distance between the two dots, and Iq is 
the zeroth-order modified Bessel function. We see that 
U12 is primarily determined by a. 

These considerations make the fit using the biquadratic 
potential straightforward, which we outline taking Si- 
MOS as an example. For the case of Fig. [TJa) the two 
dots can be regarded as basically independent so uj\ and 
uj 2 are determined from the Hubbard parameters U\ and 
U 2 using Eq. ( [lQ| ). From an appropriate U\ 2 , the inter- 
dot distance 2a is determined by Eq. (|12[ ). For the cases 
of Fig. [TJb) and (c), Eqs. (10) and |l2| provide a ini- 
tial guess from which fine tuning takes place. One key 
feature of the biquadratic confinement model, making it 
preferable for our theory, is that the minimal set of three 
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parameters entering our generalized Hubbard model (i.e. 
Ui, U2, U12) are determined by exactly three parameters 
(cji, c^2, a) of the underlying microscopic confinement 
model. Attention must be paid when transforming from 
the (/i2,Mi) to (Vr,Vl) planes: the capacitance model 
is not valid in the cases of Fig. [TJb) and (c), therefore 
we are not restricted to the restriction Eq. (m), although 
we try to choose (0.1,0.2) close to that predicted by the 
capacitance model. 

We also do the microscopic calculation with the Gaus- 
sian potential, 54 



Vo(x,y) 



-H2 exp 



exp 



(x + ai) +y 2 



J 2 



(x-a 2 ) + y 2 
J 2 



V b exp 



_ y*_ 

'l2 72 

L bx L by 



(13) 



The Gaussian potential involves seven independent pa- 
rameters. a\ and d2 denote the distance from the cen- 
ter of the quantum dot to the middle potential barrier. 
ldi and ld2 represent the size of the individual poten- 
tial well. Vb is the height of the middle potential bar- 
rier whose size is determined by 1\> X and l^y in the x- 
and ^-direction, respectively. Although simple analytical 
formulas connecting these microscopic parameters and 
the Hubbard model do not exist (particularly since the 
Hubbard model is characterized by three parameters and 
the Gaussian confinement model by seven parameters), 
the similar roles played by some parameters in the Gaus- 
sian model and the biquadratic model simplify the search. 
The parameters a\ and <22 in Eq. ( pj| should have simi- 
lar order of magnitude as a in Eq. pi), while ldi and 1^2 
for the Gaussian potential are closely related to uj\ and 
U02 in the biquadratic potential. After these parameters 
are appropriately chosen, both the height and the size 
of the central potential barrier are fine tuned to fit the 
experimental plots. 

The appropriate dielectric constant e (7.8 for Si-MOS 
and 12.375 for Si/SiGe) and the effective mass (m* = 
0.191m e ) are used throughout the microscopic calcula- 
tions. 



IV. RESULTS 
A. Si-MOS 

Here we discuss the experiments on Si-MOS devices.^ 
We start with the results from the Hubbard model, shown 
in Fig. [TJd)-(f). As discussed above, we assume that the 
capacitance model is valid in Fig. [TJa), meaning that in 
the generalized Hubbard model we need only retain pa- 
rameters Ui, U2, and Ui2- In fact, in Ref. l42l the charg- 
ing energies, namely the on-site Coulomb interactions JJ\ 
and JJii are quantitatively given as U\ = 9.8 meV and 
U2 = 11 meV. We note that dot 1 is slightly larger than 



dot 2, different from the case with symmetric double dots 
discussed in Ref. HU The gate voltages are transformed 
from the chemical potentials according to Eq. (|7|). With 
fixed U\ =9.8 meV and U2 = 11 meV, the charge stabil- 
ity diagram in this case solely depends on U12 and a fit 
to experiment can be readily obtained. We have found 
that U12 = 1.8 meV as the appropriate parameter value 
and the result is shown in Fig. [ljd) . 

To fit Fig. [jjb) and (c) we switch on the most sig- 
nificant quantum parameter in our generalized Hubbard 
model: the tunnel coupling t. We note that the val- 
ues of J7i, U2 and U12 need to be accordingly adjusted 
and we fix the ratio U1/U2 to be equal to the case of 
Fig. (TJd). The results are shown in Fig. [TJe) and (f). 
In Fig. file), U12 = 2.6 meV is intermediate compared 
to Ui = 8.36 meV and U2 = 9.38 meV which leads to 
the honeycomb shape, while a tunnel coupling of t = 0.3 
meV slightly rounds the phase boundaries near the triple 
points, in agreement with the experiment. In Fig. [TJf), 
U12 = 3.33 meV is relatively large compared to U\ = 6.53 
meV and U2 = 7.33 meV, while t = 0.4 meV. This makes 
the phase boundaries almost like diagonal parallel lines 
and the rounding effect is substantial. We note that in 
the same N block the boundaries separating different 
(711,712) states are calculated according to their respec- 
tive probabilities 45 but not the differential conductance. 
Therefore the fading effect of those lines in the experi- 
ment due to quantum fluctuations is not seen in Fig.[TJf). 
The conductance calculations involve completely differ- 
ent physics where the excited states play a role which 
is well beyond the scope of our work where we restrict 
ourselves to the ground state physics and the associated 
charge stability diagram. 

The results of the microscopic calculation with bi- 
quadratic potential [Eq. ([sj) ] are shown in Fig. [ljg)- 
(i). For Fig. [TJg) one sees excellent agreement with ex- 
periments and Hubbard model calculations. As far as 
Fig. [ljb) and (c) are concerned, the inter-dot coupling 
becomes important and Eqs. (10) and (12) cannot be 
directly applied. However, as mentioned in the above 
section they do provide a reasonable initial guess which 
refines the searching in the parameter space. Fig. [TJh) 
and (i) present the results, and good agreement with both 
the experiments and the Hubbard model calculations are 
found. We note that the coefficients (ai, 02) which trans- 
form from the (u2, Mi) to the (V#, Vl) planes are not cho- 
sen from Eq. (ITl since the capacitance model is invalid. 
Their values are provided in the caption of Fig. [I] 

We also show the results from the Gaussian potential 
[Eq. ( [l3| )] in Fig. [TJj)-(l). Very good agreement with 
their respective counterparts is evident. The Gaussian 
potential contains seven parameters which can be inde- 
pendently tuned, as such it must be noted that the data 
fitting procedure involves a large degree of freedom of 
choosing parameters which is to some extent artificial, 
especially in our case where both the precise form of the 
confinement potential and the detailed electronic configu- 
ration in the system are unknown. We have made efforts 
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FIG. 1: (Color online) Comparison of calculated charge stability diagrams and the experiments on Si-MOS system™ (a)-(c): 
experimental results reproduced from Fig. 2(a)-(c) of Ref. |3U The blue squares on the figures indicate the data range that 
we intend to fit. (d)-(f): Results from the Hubbard model, which fit panels (a), (b), and (c), respectively. Parameters: (d): 
Ui = 9.8 meV, U 2 = 11 meV, U 12 = 1.8 meV, t = 0, ai = 0.862, a 2 = 0.842. (e): Ui = 8.36 meV, U 2 = 9.38 meV, U 12 = 2.6 
meV, t = 0.3 meV, ai = 0.791, a 2 = 0.754. (f): Ui = 6.53 meV, U 2 = 7.33 meV, U 12 = 3.33 meV, t = 0.4 meV, ai = 0.710, 
a 2 = 0.638. All other parameters J e , J p , Jti, Jt 2 are zero, (g)-(i): Microscopic calculations [fitting (a)-(c)] with the biquadratic 
potential [Eq. Q]. Parameters: (g): ftcji = 0.718 meV, hw 2 = 0.905 meV, a = 52.544 nm, a x = 0.862, a 2 = 0.842. (h): 
hhj x 0.532 meV, hw 2 0.671 meV, a 37.827 nm, a± 0.772, a 2 0.732. (i): hun 0.336 meV, 0.421 meV, 

a = 28.515 nm, a\ — 0.671, a 2 = 0.608. (j)-(l): Microscopic calculations with the Gaussian potential [Eq. ( |13| )]. Parameters: 
(j) : Vb — 37 meV, a\ — 70 nm, a 2 — 56 nm, lai — 62 nm, l^ 2 — 58 nm, l bx — 15 nm, l by — 40 nm, ai — 0.885, a 2 — 0.866. (k): 
Vb — 14 meV, a\ — 36 nm, a 2 — 33 nm, Idi — 57 nm, l^ 2 — 55 nm, lb x — 18 nm, lb y — 25 nm, ol\ — 0.983, a 2 — 0.948. (1): 
Vb — 2.2 meV, a\ — 26 nm, a 2 — 20 nm, l dl — 55 nm, l d2 — 51 nm, l bx — 30 nm, l by — 55 nm, ai — 0.997, a 2 — 0.965. 
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FIG. 2: (Color online) Comparison of calculated charge stability diagrams and experiments on Si/SiGe system.^ (a)-(c): 
experimental results reproduced from Fig. 2(c)- (a) of Ref. 1431 (Note the reverse order) (d)-(f): Results from the Hubbard 
model, which fit the content enclosed by the blue frames in panels (a), (b), and (c), respectively. Parameters: (d): Ui = 8.8 
meV, U 2 6.2 meV, Ui 2 1.6 meV, t 0.09 meV, ai 0.73, a 2 0.68. (e): Ui 6.2 meV, U 2 6.1 meV, Ui 2 2.3 meV, 
t = 0.12 meV, a x = 0.653, a 2 = 0.63. (f): Ui = 4 meV, U 2 = 4 meV, U 12 = 3.2 meV, t = 0.5 meV, a x = 0.515, a 2 = 0.455. 
All other parameters J e , J p , Jti, Jt2 are zero, (g)-(i): Microscopic calculations with the biquadratic potential. Parameters: 
(g): hujx 1.397 meV, hw 2 0.723 meV, a 32.825 nm, ai 0.70, a 2 0.66. (h): fci 0.799 meV, 0.773 meV, 

a = 29.339 nm, ai = 0.675, a 2 = 0.635. (i): /kJi = 0.421 meV, ^cj 2 = 0.421 meV, a = 2.549 nm, ai = 0.515, a 2 = 0.48. (j)-(l): 
Microscopic calculations with the Gaussian potential. Parameters: (j): Vb — 10 meV, a\ — 25 nm, a 2 = 37 nm, Zdi = 30 nm, 
ld2 — 70 nm, l bx — 10 nm, lb y — 20 nm, ai — 0.95, a 2 — 0.68. (k): Vb = 9 meV, ai — 19 nm, a 2 — 18 nm, Z dl — 40 nm, l^ 2 — 43 



nm, lbx — 25 nm, 



20 nm, a\ — 0.91, a 2 — 0.86. (1): Vb — 2 meV, ai = 4 nm, a 2 = 4 nm, Z^i — 46 nm, Z^ 2 — 46 nm, 



/bcc — 5 nm, Zbj/ — 10 nm, a\ — 0.98, a 2 — 0.92. 
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to precisely reflect all reliable information that we cur- 
rently have, however these restrictions still leave several 
free degrees of freedom. In the present work we have 
shown that our theory is capable of giving a quantitative 
explanation of experiments, with all reliably known re- 
strictions incorporated. However an ultimate complete 
description of the experiment would require both im- 
provement of the experimental technique and the first- 
principles Poisson-Schrodinger approach which treats the 
exact solution to the Schrodinger equation and a real- 
istic confinement potential therein. Both of these are 
currently impractical to achieve since all the details of 
the actual potential confinement of the quantum dots are 
simply unknown for the experimental structures. 



B. Si/SiGe 

In this subsection we discuss the charge stability di- 
agram and the tunnel coupling of the Si/SiGe devices. 
In Fig. [2ja)-(c) we replicate the experimental plot from 
Fig. 2(c)- (a) of Ref. 03j The order is reversed with re- 
spect to Ref. |43] such that the tunnel coupling goes from 
weak to strong from panel (a) to (c). There are some 
noisy lines parallel to the x-axis (Vr) which slightly dis- 
locate the diamonds. They are believed to result from 
impurities in the quantum dot!^ and are not further dis- 
cussed here. The rounded boundary lines near the triple 
points mean that even for the case of Fig. [2|a) one needs 
to take the tunnel coupling into account. Moreover, the 
asymmetry between the two dots is most pronounced in 
Fig. |2|a) but is substantially suppressed in Fig. [2^b) and 
(c), as can be seen directly by reflecting the figure with 
respect to the diagonal line Vl~Vr — const.. We suspect 
that it is due to the fundamental difference in the way the 
Si-MOS and Si/SiGe devices are fabricated, therefore un- 
like the Si-MOS case where we fix the ratio U1/U2, in the 
Si/SiGe case we use independent U\ and U2 values that 
are appropriate for the given charge stability diagrams. 

Fig. [2jd)-(f) present the results from the Hubbard 
model. In Fig.^Jd), the asymmetry is evident as V\ = 8.8 
meV and U2 = 6.2 meV. A finite t = 0.09 meV makes the 
boundary lines slightly rounded. In Fig. |2|e) the values 
of U\ and U2 are already quite similar (U\ = 6.2 meV, 
U2 = 6.1 meV), implying a rather symmetric honeycomb 
pattern, t — 0.12 meV, which is larger than that of 
Fig. [2|d), makes the boundary line even more rounded. 
As stated in Ref. [43j in Fig. [2^f ) the double dot system 
effectively becomes a single, large quantum dot and our 
fitted parameters U\ = U2 = 4 meV, and rather large 
values of /7i2 = 3.2 meV and t = 0.5 meV are consistent 
with this argument. We note that in Fig. |5|f) the (0, 2) 
and (2, 0) components of the stability diagram appear 
out of the range shown in the figure. This can be un- 
derstood since, as the two dots effectively merge due to 
very low central barrier, the (1,1) state dominates most 
of the diagram and one must use substantially different 
chemical potentials on the two dots in order to achieve 



double occupancy on one but not the other of the two 
dots. 

We show the results of the microscopic calculation from 
the biquadratic model potential and the Gaussian poten- 
tial in Fig. [2jg)-(i) and Fig. |2jj)-(l), respectively. The 
overall features are observed to be very similar as the pa- 
rameters are tuned, but there are differences in the fine 
details. For example in Fig. [2|j) the boundaries between 
different (774,77,2) states within the same N block [e.g. 
(1,1) and (0,2)] are tilted, and in Fig. ^i) and (1) the 
spacing of the diagonal "parallel lines" are observed to 
be quite different from the experiments. This indicates 
that the model confinement potential is sometimes insuf- 
ficient to mimic the true potential that electrons would 
feel in the quantum dot system. To resolve this problem 
one must calculate the electrostatic confinement poten- 
tial according to the design of the devices using, for ex- 
ample, the Poisson equation. Since obtaining the precise 
connection between gates and the electron gas localized 
inside the quantum dots is difficult, one would not ex- 
pect that this problem can be straightforwardly resolved 
along this line. Further study of this issue is needed. 

The tunnel coupling t is crucial in the double quan- 
tum dot system since it implies the exchange interaction, 
which is essential for many proposals of spin-based quan- 
tum computationP^Sl While the exchange interaction 
can be directly measured experimentally from the sin- 
glet/triplet level splitting, 42 the tunnel coupling t is of 
fundamental theoretical importance as it appears in the 
Hamiltonian matrix. As discussed in Ref. 45, there are 
at least two indirect ways to measure t: one is through 
the curvature of the boundary lines, as discussed above 
in the context of stability diagrams; the other is accord- 
ing to the width of the probability crossover near the 
phase boundaries within a subspace of the charge stabil- 
ity diagram, for example the (1,1) and (0, 2) probability 
crossing. The latter strategy has been used in Refs. l43l 
and [62] the crossover of (774 + 1,772) and (771,772 + 1) 
states leads to an obvious change in the tunneling current 
which is directly measured by the quantum point contact. 
Then employing a two-level model (which is essentially 
a subspace of our more general model) one extracts the 
t values. Fig. 4 of Ref. [43] shows such a crossover of the 
occupancy and its inset shows the extracted values of t as 
a function of the middle gate voltage Vm, whose fitting 
line is reproduced as the black solid line in the inset of 
Fig. [3] As the absolute value of Vm increases, the central 
potential barrier is raised and t decreases approximately 
exponentially. 

Here we present a microscopic calculation which shows 
similar behavior and compare that to the experiments. 
We use the biquadratic model potential as it contains 
fewer free parameters and the inter-dot distance has 
a simple relation to the height of the potential bar- 
rier and the level spacing of the harmonic oscillator: 
2a = 2^2ht)/m/uJo. Since the precise microscopic pa- 
rameters cannot be known from the experiments, we use 
a range of different harmonic oscillator level spacings, 
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FIG. 3: (Color online) Tunnel coupling t calculated from the biquadratic model potential as functions of the effective central 
gate voltage Vb = hb/(—e), plotted on the semi-log scale. Here hb is the height of the central potential barrier. Since t varies 
over the whole range of the charge stability diagram, we choose the chemical potentials at the center of (1,1) phase. As the 
absolute value of Vb decreases, t increases approximately exponentially. The harmonic oscillator energy level spacings of the 
single-dot confinements are fixed as ftwo = 1.319meV (red circles), huuo = 0.989meV (blue squares), and ftwo = 0.659meV (green 
triangles). The inter-dot distance 2a = 2y/2hb/m/ujQ. Inset: Overlaying comparisons between t extracted from experiment as 
a function of middle gate voltage Vm (black solid line, bottom x-scale reproduced from Fig. 4 inset of Ref. 43) and t calculated 
from the microscopic model as a function of effective central gate voltage Vb ( huoo = 0.989meV, blue squares, top x-scale). The 
y-axis is shared by the two but note that the x-scales differ by a proportionality constant and a shift. 



assuming the two dots to be symmetric, and study the 
tunnel coupling t as a function of the height of the central 
potential barrier h^. The results are shown in the main 
panel of Fig. [3j the results of three values of harmonic os- 
cillator level spacing are shown as different symbols with 
hwo = 1.319meV (red circles), Hujq = 0.989meV (blue 
squares), and Huq = 0.659meV (green triangles). Since 
the calculated value of t changes over the entire charge 
stability diagram, we focus on the center point of the 
(1,1) component of the charge stability diagram and com- 
pare all cases at that point. To facilitate the comparison 
to experiments we convert the central barrier height to 
the effective gate voltage Vb = hb/(— e). On the semi-log 
scale the data points show an approximately exponential 
behavior. The exponential behavior can again be under- 
stood from the local nature of the electron wave function: 
the overlap between two well localized wave function at a 
distance decreases rapidly as they become either further 
apart or more localized. 

We compare our results calculated at huo = 0.989 meV 
(blue squares) to experiments in the inset of Fig. [3] Dif- 
ficulties arise since, although the experimental middle 
gate voltage Vm plays a similar role as the central po- 
tential barrier of the model potential hb in our calcu- 
lation, the quantitative connection between the realistic 
confinement potential and our model potential is in gen- 
eral unknown. First, the realistic confinement potential 
depends on details of the fabrication and the geometry 
of the specific device, which is difficult to quantitatively 
describe. Second, the multi-electron multi-band feature 



in the real situation further complicates the problem. 
Therefore, in our situation, we shall only assume that Vm 
and Vb [= hb/(—e)] are related linearly without quanti- 
fying the linear coefficients. In the inset of Fig. [3] the 
theoretical results are shown as points, utilizing the top 
x-scale, while the experimental fitting line is shown as 
the solid line, using the bottom x-scale. The y-scale is 
shared by both and the x-scales differ by a proportion- 
ality constant and a shift as discussed above. Caution 
must be taken when interpreting this figure: the perfect 
fit does not mean that the microscopic parameters used 
here are the correct ones for the corresponding experi- 
ment because of our many simplifying approximations. 
The purpose is only to demonstrate that our method is 
capable of producing the exponential behavior, assuming 
that Vm and Vb are related linearly. The precise connec- 
tion between them is still uncertain. In fact each group 
of the three hwo points can be appropriately rescaled to 
fit the experimental curve. Therefore we only argue that 
our theory produces the exponential behavior along with 
appropriate, albeit rough, simplifications, and leave the 
precise first-principles treatment for future studies. Such 
a study, which is well beyond the scope of our minimal 
model, will not only require intensive computational ef- 
forts but also necessitate much better characterization of 
the experimental dots than is currently available. 
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V. CONCLUSION 

In conclusion, we have applied the generalized two-site 
Hubbard model approach previously introduced in the 
GaAs symmetric double dot system to the Si asymmetric 
double dot system. We have calculated the charge sta- 
bility diagrams from the generalized Hubbard model and 
the microscopic theory with two different models of con- 
finement potentials (biquadratic and Gaussian forms), 
and have found qualitative and quantitative agreement 
between them and with the experiments on Si-MOS and 
Si/SiGe devices. We have explained in detail the simplifi- 
cations that we have made to apply our theory, which in- 
volves the reduction to the effective two-electron regime^ 
and approximation to the gate voltages. We have shown 
that our proposed theoretical Hubbard model approach, 
along with appropriate approximations, provides state- 
of-the-art quantitative description to the silicon double 
quantum dot system. In addition, we have studied in par- 
ticular the tunnel coupling using the microscopic calcu- 
lation with biquadratic model potential and have shown 
that the experimentally observed exponentially decaying 
behavior with the central potential barrier height is con- 
sistent with our theory. 

The qualitative success of our minimal theoretical 
model should motivate further investigations. First, the 
full understanding of the multi-electron nature of silicon 
quantum dot devices requires a multi-band calculation.^ 
In our work we have qualitatively argued on physical 
grounds that the multi-electron effect imposes simply a 
renormalization of Hubbard parameters in the general- 
ized one-band effective Hubbard model, 49 but a quanti- 
tative assessment is lacking. It is important to under- 
stand in-depth how the Hubbard parameters change as 
one translates between different multi-electron regimes, 
in particular those corresponding to quantum fluctua- 
tions which may have non-trivial dependence on the con- 
figuration of the base Fermi sea. Second, a first-principles 
Poisson-Schrodinger theory would be helpful in deter- 



mining the microscopic confinement potential. The de- 
ficiency of current experiments is that the confinement 
potential is subject to fabrication-induced disorder and 
it is hard to extract the precise form of confinement 
potential reliably. Beyond the Hund-Mulliken approx- 
imation, one can in principle obtain exact solutions to 
the Schrodinger equation for any arbitrary confinement 
potential, although for some calculations the cost is ex- 
pensive. However a complete quantitative understanding 
will be impossible unless one acquires knowledge beyond 
the model potential and the effective two-electron regime. 
Some of the above problems are intractable for current 
technologies. However, these warrant further study. A 
key point to appreciate in this context is the fact that no 
two sets of double dots are identical from the perspective 
of microscopic details even when they are fabricated from 
the same underlying two-dimensional electron structures 
using the same nominal lithographic protocols. This is 
due to many reasons, but a primary one is the invariable 
presence of unintentional (and therefore, unknown) back- 
ground disorder in the host semiconductor, particularly 
at the interface where the dots reside. It may there- 
fore be useless trying to go much beyond our effective 
minimal model at this stage of materials development in 
solid state spin qubits since no generic quantitative the- 
ory is feasible until the characterization of the quantum 
dots and their underlying materials science become much 
more advanced and reproducible than they currently are. 
We believe that our minimal model is perhaps the most 
practical quantum generalization of the classical capaci- 
tance model that is feasible at this stage of development 
of solid state spin qubits in semiconductors. 

Acknowledgements 

We thank M. A. Eriksson, J. P. Kestner and Q. Li 
for helpful discussions. This work is supported by LPS- 
CMTC. 



1 D. Loss and D. P. DiVincenzo, Phys. Rev. A 57, 120 
(1998). 

2 D. P. DiVincenzo, D. Bacon, J. Kempe, G. Burkard, and 
K. B. Whaley, Nature (London) 408, 339 (2000). 

3 J. Levy, Phys. Rev. Lett. 89, 147902 (2002). 

4 F. H. L. Koppens, C. Buizert, K. J. Tielrooij, I. T. Vink, 
K. C. Nowack, T. Meunier, L. P. Kouwenhoven, and L. M. 
K. Vandersypen, Nature (London) 442, 766 (2006). 

5 B. E. Kane, Nature (London) 393, 133 (1998). 

6 R. Vrijen, E. Yablonovitch, K. Wang, H. W. Jiang, A. 
Balandin, V. Roychowdhury, T. Mor, and D. DiVincenzo, 
Phys. Rev. A 62, 012306 (2000) 

7 M. Friesen, P. Rugheimer, D. E. Savage, M. G. Lagally, D. 
W. van der Weide, R. Joynt, and M. A. Eriksson, Phys. 
Rev. B 67, 121301 (R) (2003) 

8 J. M. Taylor, H.-A. Engel, W. Diir, A. Yacoby, C. M. Mar- 



cus, P. Zoller, and M. D. Lukin, Nat. Phys. 1, 177 (2005). 
9 J. M. Kikkawa and D. D. Awschalom, Phys. Rev. Lett. 80, 
4313 (1998). 

10 S. Amasha, K. MacLean, Iuliana P. Radu, D. M. Zumbiihl, 
M. A. Kastner, M. P. Hanson, and A. C. Gossard, Phys. 
Rev. Lett. 100, 046803 (2008) 

11 F. H. L. Koppens, K. C. Nowack, and L. M. K. Vander- 
sypen, Phys. Rev. Lett. 100, 236802 (2008) 

12 C. Barthel, J. Medford, C. M. Marcus, M. P. Hanson, and 
A. C. Gossard, Phys. Rev. Lett. 105, 266808 (2010) 

13 H. Bluhm, S. Foletti, I. Neder, M. Rudner, D. Mahalu, V. 
Umansky, and A. Yacoby, Nat. Phys. 7, 109 (2011). 

14 A. M. Tyryshkin, S. A. Lyon, A. V. Astashkin, and A. M. 
Raitsimring, Phys. Rev. B 68, 193207 (2003). 

15 A. Morello, J. J. Pla, F. A. Zwanenburg, K. W. Chan, K. 
Y. Tan, H. Huebl, M. Mottonen, C. D. Nugroho, C. Yang, 



11 



J. A. van Donkelaar, A. D. C. Alves, D. N. Jamieson, C. 
C. Escott, L. C. L. Hollenberg, R. G. Clark, and A. S. 
Dzurak, Nature (London) 467, 687 (2010). 

16 C. B. Simmons, J. R. Prance, B. J. Van Bael, T. S. Koh, 
Z. Shi, D. E. Savage, M. G. Lagally, R. Joynt, M. Friesen, 
S. N. Coppersmith, and M. A. Eriksson, Phys. Rev. Lett. 
106, 156804 (2011). 

17 M. Xiao, M. G. House, and H. W. Jiang, Phys. Rev. Lett. 
104, 096801 (2010) 

18 J. R. Petta, A. C. Johnson, J. M. Taylor, E. A. Laird, A. 
Yacoby, M. D. Lukin, C. M. Marcus, M. P. Hanson, A. C. 
Gossard, Science 309, 2180 (2005). 

19 S. Foletti, H. Bluhm, D. Mahalu, V. Umansky, and A. 
Yacoby, Nat. Phys. 5, 903 (2009) 

20 E. A. Laird, J. M. Taylor, D. P. DiVincenzo, C. M. Marcus, 
M. P. Hanson, and A. C. Gossard, Phys. Rev. B 82, 075403 
(2010). 

21 I. van Weperen, B. D. Armstrong, E. A. Laird, J. Medford, 

C. M. Marcus, M. P. Hanson, and A. C. Gossard, e-print 
|arXiv:1101.5654| (unpublished) . 

22 C. Tahan, M. Friesen, and R. Joynt, Phys. Rev. B 66, 
035314 (2002) 

23 R. de Sousa and S. Das Sarma, Phys. Rev. B 68, 115322 
(2003) 

24 L. Fedichkin and A. Fedorov, Phys. Rev. A 69, 032311 
(2004) 

25 I. Appelbaum, B. Huang, and D. J. Monsma, Nature (Lon- 
don) 447, 295 (2007) 

26 E. P. Nordberg, G. A. Ten Eyck, H. L. Stalford, R. P. 
Muller, R. W. Young, K. Eng, L. A. Tracy, K. D. Childs, 
J. R. Wendt, R. K. Grubbs, J. Stevens, M. P. Lilly, M. 
A. Eriksson, and M. S. Carroll, Phys. Rev. B 80, 115331 
(2009) 

27 S. Goswami, K. A. Slinker, M. Friesen, L. M. McGuire, J. 
L. Truitt, C. Tahan, L. J. Klein, J. O. Chu, P. M. Mooney, 

D. W. van der Weide, R. Joynt, S. N. Coppersmith, and 
M. A. Eriksson, Nat. Phys. 3, 41 (2007) 

28 M. G. Borselli, R. S. Ross, A. A. Kiselev, E. T. Croke, K. 
S. Holabird, P. W. Deelman, L. D. Warren, I. Alvarado- 
Rodriguez, I. Milosavljevic, F. C. Ku, W. S. Wong, A. E. 
Schmitz, M. Sokolich, M. F. Gyure, and A. T. Hunter, 
Appl. Phys. Lett. 98, 123118 (2011) 

29 S. Das Sarma, R. de Sousa, X. Hu, and B. Koiller, Sol. St. 
Comm. 133, 737 (2005) 

30 D. Culcer, L. Cywihski, Q. Li, X. Hu, and S. Das Sarma, 
Phys. Rev. B 80, 205302 (2009); 82, 155312 (2010) 

31 Q. Li, L. Cywihski, D. Culcer, X. Hu, and S. Das Sarma, 
Phys. Rev. B 81, 085313 (2010). 

32 B. Hu and C. H. Yang, Phys. Rev. B 80, 075310 (2009) 

33 W. H. Lim, F. A. Zwanenburg, H. Huebl, M. Mottonen, 
K. W. Chan, A. Morello, and A. S. Dzurak, Appl. Phys. 
Lett. 95, 242102 (2009) 

34 M. Xiao, M. G. House, and H. W. Jiang, Appl. Phys. Lett. 
97, 032103 (2010) 

35 L. J. Klein, D. E. Savage, and M. A. Eriksson, Appl. Phys. 
Lett. 90, 033103 (2007) 

36 C. B. Simmons, M. Thalakulam, N. Shaji, L. J. Klein, 
H. Qin, R. H. Blick, D. E. Savage, M. G. Lagally, S. N. 
Coppersmith, and M. A. Eriksson, Appl. Phys. Lett. 91, 
213103 (2007) 

37 N. Shaji, C. B. Simmons, M. Thalakulam, L. J. Klein, H. 
Qin, H. Luo, D. E. Savage, M. G. Lagally, A. J. Rimberg, 
R. Joynt, M. Friesen, R. H. Blick, S. N. Coppersmith, and 
M. A. Eriksson, Nat. Phys. 4, 540 (2008) 



38 C. B. Simmons, T. S. Koh, N. Shaji, M. Thalakulam, L. 
J. Klein, H. Qin, H. Luo, D. E. Savage, M. G. Lagally, 
A. J. Rimberg, R. Joynt, R. Blick, M. Friesen, S. N. Cop- 
persmith, and M. A. Eriksson, Phys. Rev. B 82, 245312 
(2010) 

39 H. W. Liu, T. Fujisawa, Y. Ono, H. Inokawa, A. Fujiwara, 
K. Takashina, and Y. Hirayama, Phys. Rev. B 77, 073310 
(2008) 

40 M. Pierre, R. Wacquez, B. Roche, X. Jehl, M. Sanquer, 
M. Vinet, E. Prati, M. Belli, and M. Fanciulli, Appl. Phys. 
Lett. 95, 242107 (2009) 

41 C. Livermore, C. H. Crouch, R. M. Westervelt, K. L. 
Campman, and A. C. Gossard, Science 274, 1332 (1996). 

42 N. S. Lai, W. H. Lim, C. H. Yang, F. A. Zwanenburg, 
W. A. Coish, F. Qassemi, A. Morello, and A. S. Dzurak, 
e-print |arXiv:1012.1410| (unpublished) . 

43 C. B. Simmons, M. Thalakulam, B. M. Rosemeyer, B. J. 
Van Bael, E. K. Sackmann, D. E. Savage, M. G. Lagally, R. 
Joynt, M. Friesen, S. N. Coppersmith, and M. A. Eriksson, 
Nano Lett. 9 3234 (2009). 

44 L. DiCarlo, H. J. Lynch, A. C. Johnson, L. I. Childress, K. 
Crockett, C. M. Marcus, M. P. Hanson, and A. C. Gossard, 
Phys. Rev. Lett. 92, 226801 (2004). 

45 S. Yang, X. Wang, and S. Das Sarma, Phys. Rev. B 83, 
161301(R) (2011). 

46 W. G. Wiel, S. De Franceschi, J. M. Elzerman, T. Fujisawa, 
S. Tarucha, and L. P. Kouwenhoven, Rev. Mod. Phys. 75, 
1 (2003). 

47 R. Hanson, L. P. Kouwenhoven, J. R. Petta, S. Tarucha, 
and L. M. K. Vandersypen, Rev. Mod. Phys. 79, 1217 
(2007). 

48 J. M. Elzerman, R. Hanson, J. S. Greidanus, L. H. Willems 
van Beveren, S. De Franceschi, L. M. K. Vandersypen, 
S. Tarucha, and L. P. Kouwenhoven, Phys. Rev. B 67, 
161308(R) (2003). 

49 X. Wang, S. Yang, and S. Das Sarma, e-print 
larXiv:1104.549Tl 

50 D. Culcer, X. Hu, and S. Das Sarma, Phys. Rev. B 82, 
205315 (2010). 

51 L. Gaudreau, S. A. Studenikin, A. S. Sachrajda, P. Za- 
wadzki, A. Kam, J. Lapointe, M. Korkusinski, and P. 
Hawrylak, Phys. Rev. Lett. 97, 036807 (2006). 

52 M. Korkusinski, I. P. Gimenez, P. Hawrylak, L. Gaudreau, 
S. A. Studenikin, and A. S. Sachrajda, Phys. Rev. B 75, 
115301 (2007). 

53 G. Burkard, D. Loss, and D. P. DiVincenzo, Phys. Rev. B 
59, 2070 (1999). 

54 X. Hu and S. Das Sarma, Phys. Rev. A 61, 062301 (2000). 

55 X. Hu and S. Das Sarma, Phys. Rev. A 64, 042312 (2001). 

56 R. de Sousa, X. Hu, and S. Das Sarma, Phys. Rev. A 64, 
042307 (2001). 

57 I. Puerto Gimenez, M. Korkusinski, and P. Hawrylak, 
Phys. Rev. B 76, 075336 (2007). 

58 E. Nielsen, R. W. Young, R. P. Muller, and M. S. Carroll, 
Phys. Rev. B 82, 075319 (2010) 

59 M. Helle, A. Harju, and R. M. Nieminen, Phys. Rev. B 72, 
205329 (2005) 

60 J. Pedersen, C. Flindt, N. A. Mortensen, and A.-P. Jauho, 
Phys. Rev. B 76, 125323 (2007) 

61 N. T. T. Nguyen and S. Das Sarma, Phys. Rev. B (in 
press) , |arXiv: 1103.0767| 

62 T. Hatano, M. Stopa, and S. Tarucha, Science 309, 268 
(2005). 



